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We present a large- N theory for the differential conductance, dl/dV, in Kondo systems measured 
via scanning tunneling spectroscopy. We demonstrate that quantum interference between tunneling 
processes into the conduction band and into the magnetic /-electron states is crucial in determining 
the experimental Fano lineshape of dl/dV. This allows one to uniquely extract the Kondo coupling 
and the ratio of the tunneling amplitudes from the experimental dl/dV curve. Finally, we show that 
dl/dV directly reflects the strength of the antiferromagnetic interaction in Kondo lattice systems. 
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Recent progress in scanning tunneling spectroscopy 
(STS) techniques have made it possible for the first 
time to measure the differential conductance, dl/dV, in 
heavy- fermion compounds [lj. These materials, whose 
essential ingredients are a (Kondo) lattice of magnetic 
moments that is coupled to a conduction band [2j, ex- 
hibits a variety of puzzling phenomena, ranging from 
non-Fermi-liquid behavior to unconventional supercon- 
ductivity [3|- Their microscopic origin likely lies in the 
competition between an antiferromagnetic ordering of 
the magnetic moments, and their screening by conduc- 
tion electrons [2], though no theoretical consensus has 
emerged as yet [4j|. STS experiments, by providing in- 
sight into the local electronic structure [1] of heavy- 
fermion materials, might hold the key to understand- 
ing their complex properties. The theoretical challenge 
in the interpretation of the differential conductance in 
Kondo lattice systems [5], and around single Kondo im- 
purities arises from the quantum interference be- 
tween electrons tunneling into the conduction band and 
into the magnetic /-electron states. While dl/dV for a 
single Kondo impurity has been successfully described 
Q using a phenomenological form derived by Fano [TlJ , a 
microscopic understanding of how the interplay between 
the strength of the Kondo coupling, the interaction be- 
tween the magnetic moments, the electronic structure of 
the screening conduction band, and quantum interference 
determines the dl/dV lineshape, is still lacking. 

In this Letter, we address this issue within the frame- 
work of a large- N theory and identify the microscopic 
origin for the form of the differential conductance not 
only around single Kondo impurities but also in Kondo 
lattice systems. In particular, we demonstrate that the 
lineshape as well as the spatial dependence of dl/dV sen- 
sitively depend on the particle-hole asymmetry of the 
(screening) conduction band, as well as the quantum in- 
terference between the two tunneling paths. For a single 
Kondo impurity, this sensitivity allows one to uniquely 
extract the Kondo coupling, J, as well as ratio of the 
tunneling amplitudes into the conduction band and mag- 
netic /-electron states, t c and tf, respectively, from the 
experimental STS data. In addition, for a Kondo lattice, 



the dl/dV lineshape provides insight into the strength 
of the interaction between the magnetic moments. Due 
to quantum interference effects, which can lead to a re- 
versal in the asymmetry of the dl/dV lineshape already 
for small changes in tf/t c , the differential conductance 
is in general qualitatively different from the local density 
of states (LDOS) of either the conduction band or the 
/-electron states. However, once the pertinent parame- 
ters are extracted from a theoretical fit, we can predict 
the frequency and spatial dependence of the LDOS for 
both bands, as well as the electronic correlations between 
them, thus providing important insight into the complex 
electronic structure of Kondo systems. 
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FIG. 1: (color online) Tunneling paths from the STS tip into 
the conduction electron and /-electron states with tunneling 
amplitudes t c and tf, respectively. 

We start by considering the differential conductance in 
a system with a single Kondo impurity, whose Hamilto- 
nian is given by 

H = - *rr'4 )a C r> + JSg-S^, (1) 

r,r',cr 

where t rr / is the fermionic hopping element between sites 
r and r' of the conduction band, c\ a (c r>(T ) creates (an- 
nihilates) a fermion with spin a at site r, and the sums 
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run over all sites of the conduction band. and 
are the spin operators of the Kondo impurity and the 
conduction electrons at site R, respectively, and J > 
is the Kondo coupling. To describe the Kondo screening 
of the magnetic impurity, we employ a large-N expansion 
12 lH in which is generalized to SU(N) and repre- 
sented via Abrikosov pseudofermions , / m • These obey 
the constraint J^ m =i n fmfm = 1 where N = 2S + 1 is 
the spin degeneracy of the magnetic impurity. This con- 
straint is enforced by means of a Lagrange multiplier ef, 
while the exchange interaction in Eq.(pQ) is decoupled via 
a hybridization field, s. On the saddle point level, Sf and 
s 2 are obtained by minimizing the effective action [13j . 
Finally, the tunneling process into a conduction electron 
state at r and the /-electron state at R, as schematically 
shown in Fig. [TJ is described by 



U T = ^2t c cl a d a + tff^da + H.c. , 



(2) 



where d a destroys a fermion in the STS tip. The total 
current flowing from the tip into the systems is (l9( 



n Re J ^ 



t c G«(u) + t f G? 3 (") -(3) 



I(V) : 

Here 

G K (w) = [i - S'MtrVoMIi - tfl°(w)] _1 (4) 

is the full Keldysh Greens function matrix, t is the sym- 
metric hopping matrix with £ 12 = t c , £13 = £/, and zero 
otherwise. g r (oj) and fo(oj) are the retarded and Keldysh 
Greens function matrices of the Kondo system with 

fo(uj) = 2i(l-2h F (cu))lm[g r (u)} ; 

(9t(u) \ 

g r (u) = g r Jrw) g r cf(r,-R,Lj) , (5) 

V g r fc (H,r,uj) g T ff(R, R, w)) 

where g r t ,g r cc and g r ^ are the local Greens functions of 
the tip, conduction and /-electron states, respectively, 
and # /c (R,r,r) = -(T r /^(r)c r (0)). Moreover, 
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with (rip) being the Fermi-distribution function of the 
tip (/- and c-electron states). While the results shown 
below are obtained from Eq.© via differentiation, it is 
instructive to consider the leading order contributions to 
dl/dV in the weak-tunneling limit, t c: tf —> 0, given by 



dI(V) - 2 -^N t [tlN c (r,V)+t}N f (K,V) 



dV 



2t c t f N cf (r,R, V)} 



(7) 



with Nt, N c and Nf being the density of states of the tip, 
conduction and /-electron states, respectively, and N c f = 
—Img^/ir. The last term in Eq.([7]) describes quantum 
interference processes between the two tunneling paths, 
which, as we show below, are crucial in determining the 
lineshape of the differential conductance. 
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FIG. 2: (color online) (a) Experimental dl/dV curve of 
Ref. at the site of a Co atom on a Au(lll) surface together 
with a theoretical fit from Eq.® with N = 4, t f /t c = 0.066, 
t c = 1 meV, J — 1.39 eV, s 250 meV, e f = 19 meV and 
N t = 1/eV. A constant background was subtracted from the 
experimental data, (b) conduction electron LDOS, N c (uo) (c) 
/-electron LDOS N f (u), and (d) N cf (u) at R. (e) dl/dV at 
a distance of r = 3ao from the Co atom for tf = 0. Parame- 
ters in (b)-(e) are the same as in (a), (f) Experimental dl/dV 
curve of Ref. 0] at r — 5 A from the Co atom. 

In Fig. [2fa) we present the experimental dl/dV data 
of Ref. [9[ for a tip positioned above a Co atom on a 
Au(lll) surface together with a theoretical fit obtained 
from Eq.(|3|). Here, tunneling into the conduction band 
involves only the state at R, i.e., r = R. The peak and 
dip in dl/dV are a direct signature of the hybridization 
between the conduction band and the /-electron state 
of the Kondo impurity and are commonly referred to as 
the Kondo resonance. As input parameters, we took the 
screening conduction band to be given by the Au(lll) 
surface states possessing a triangular lattice structure 
with t = 1.3 eV and fi = -7.34 eV [20], and used N = 4 
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to describe the S = 3/2 Co spin. The theoretical dl/dV 
curve is then solely determined by J and tf/t c , which 
control the width of the dip and the overall asymme- 
try of dl/dV, respectively. By performing an extensive 
survey, we found that there exists a unique set of pa- 
rameters, J = 1.39 eV and tf/t c = 0.0066, that yield 
the good quantitative agreement between the theoretical 
and experimental data shown in Fig. [2(a). We note that 
while the STS tip is positioned above the Co atom, tf/t c 
is small, likely reflecting the suppression of the tunnel- 
ing process into the /-electron state by Coulomb effects. 
Moreover, once J is obtained from the fit, we can com- 
pute the LDOS of the conduction and /-electron states, 
which are presented in Figs. [2(b) and (c), respectively, as 
well as the electronic correlations between the two bands, 
as reflected by N c f shown in Fig. [2(d). The lineshape of 
either LDOS (or of their superposition) is qualitatively 
different from that of dl/dV, demonstrating the impor- 
tance of quantum interference in determining the latter. 
Finally, as the STS tip is moved away from the Co atom, 
direct tunneling into the /-electron state becomes sup- 
pressed and tf — »• 0. Therefore, in Fig. [2(e), we present 
the theoretical dl/dV curve with tf = at a distance of 
r = 3ao from the Co atom. We note that while tf = 
the asymmetry of dl/dV is now the same as that at the 
site of the Co atom, and qualitatively agrees with the ex- 
perimental dl/dV curve at r = bA [6] shown in Fig. [2(f), 
demonstrating the consistency of our approach. A more 
quantitative fit would require an extensive spatial survey 
of dl/dV away from the Co atom. 

The asymmetry of the dl/dV lineshape is determined 
by two microscopic properties: the particle- hole asym- 
metry of the screening conduction band, and the ratio 
of the tunneling amplitudes, tf/t c . To demonstrate this 
dependence, we present in Fig. [3] the evolution of dl/dV 
with increasing ratio tf/t c . To contrast and complement 
the results shown in Fig.[2j we take N = 2, corresponding 
to a spin- 1/2 moment, and consider a conduction band 
on a square lattice with t = 0.5Eo and \i — — 1.809-Erj. 
The resulting circular Fermi surface with Fermi wave- 
length \p — 10ao is representative of the Au(lll) and 
Cu(lll) surfaces states employed in Refs.jil. lil [To|. For 
tf = [solid line in Fig. [3(a)], dl/dV exhibits a Kondo 
resonance whose asymmetry is opposite to the experi- 
mentally observed one shown in Fig. H(a). The asym- 
metry of dl/dV for tf = is a direct consequence of 
the particle-hole asymmetry of the conduction band. In- 
deed, reversing the latter via /i — ^ — /i, also leads to a 
reversal of the asymmetry in dl/dV, as shown by the 
dashed line in Fig. [3(a). Moreover, with increasing tf/t c , 
the height of the peak on the negative energy side, as well 
as the width of the dip in dl/dV decrease while its min- 
imum shifts to lower energies [see Fig. [3(b)], leading to 
an almost symmetric dl/dV curve for tf/t c = 0.062 [see 
Fig. [3(c)]. Increasing tf/t c even further [see Fig. [3(d)] 
now reverses the asymmetry in dl/dV, yielding a peak on 
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FIG. 3: (a) - (f) dl/dV at r = R as a function of energy 
for N = 2, J = 0.5£ , N t 1.0/E , t c 0.001£ , s f 
0.00520^o, s = 0.0847^o and different values of t f /t c . Dashed 
line in (a) represents dl/dV for a conduction band with a 
reversed particle-hole (rph) asymmetry. 



the positive energy side, and a minimum at slightly neg- 
ative energies. The asymmetry of the dl/dV lineshape is 
now similar to that observed experimentally. However, in 
contrast to the case of a spin-3/2 moment (N = 4) con- 
sidered in Fig. [2(a), the minimum in dl/dV is located 
at negative energies for a spin- 1/2 moment (N = 2). 
Indeed, for N = 2, no fit to the experimental data of 
Fig. [21(a) can be obtained. This demonstrates that the 
differential conductance directly reflects the spin of the 
screened magnetic moment. Increasing tf/t c even fur- 
ther [see Fig. [3(f)] leads to an increase in the height of 
the peak and a widening of the dip. 

We next turn to the discussion of the differential con- 
ductance in a Kondo lattice system, whose complex prop- 
erties are determined by the competition between the 
Kondo screening of the magnetic moments and their 
antiferromagnetic ordering [2]. Its Hamiltonian is ob- 
tained by appropriately extending Eq.(pQ), and adding 
the term Hi = J2 r r' ^m-'S^S^ representing the anti- 
ferromagnetic interaction between the moments. Here, 
we take I r y > to be non-zero for nearest-neighbor sites 
only. Using again an Abrikosov pseudo-fermion represen- 
tation of , the Hamiltonian is decoupled by introduc- 
ing the spatially uniform mean-fields [17| s = J(/J a c r,a) 
and xo = I(fr afr^a), where the latter is a measure 
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FIG. 4: Evolution of dl/dV in a Kondo lattice with tf/t c 
for N = 2, J = 0.5E , N t = 1/E , t c = 0.001£ , and (a) 
- ( c ) I/J = 0.001 with s f = 0.0012^o, s = 0.0485£ , Xo = 
0.00017^o, and (d) - (f) I/J = 0.015 with e f = 0.00094£ , 
s = 0.0480^o and xo = 0.00259£ . 



of the magnetic correlations in the system. The con- 
straint (rif) = 1 is enforced via a local on-site energy, 

Xlr £ ffr,afr,a- 

The magnetic interactions in the (screened) Kondo lat- 
tice have a profound effect on the form of the differen- 
tial conductance, as shown in Fig. |4] where we present 
the evolution of dl/dV with tf/t c for two different mag- 
netic interaction strengths, I/J = 0.001 (left column) 
and I/J = 0.015 (right column). While for I/J = 0.001, 
dl/dV exhibits a hard hybridization gap for all values of 
tf/t c , only a suppression of the differential conductance 
around the Fermi energy is found for I / J = 0.015. How- 
ever, in both cases, dl/dV exhibits a peak on the negative 
energy side (indicated by arrows), which arises from the 
Van Hove singularity of the large (hybridized) Fermi sur- 
face. This peak is first suppressed with increasing tf/t c 
[see Figs.H^b) and (e)], but then reemerges together with 
a second peak [indicated by arrows in Figs.|4jc) and (f)], 
which is the precursor of the emerging /-electron band. 
This second peak is centered around the Fermi energy 
for I/J = 0.001, but is located at positive energies for 
I/J = 0.015. In the latter case, we also find a shift of 
the minimum in dl/dV to lower energies with increasing 
tf/tc- This qualitative difference in the differential con- 
ductance thus provides direct insight into the strength of 



the antiferromagnetic interactions. 

In summary, we have presented a large- TV theory for 
the differential conductance in Kondo systems. We 
demonstrated that quantum interference between tun- 
neling paths is crucial in explaining the experimentally 
observed Fano lineshape of dl/dV. This allows one to 
uniquely extract the Kondo coupling as well as the ra- 
tio of the tunneling amplitudes from the experimental 
dl/dV curve. Finally, we showed that dl/dV reflects the 
strength of the antiferromagnetic interaction in Kondo 
lattice systems. 
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